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Abstract. We explain how the constraints from present experimental data can 
be used to obtain the nPDF in the framework of LO DGLAP evolution. We 
will also compare the only two available sets of this type and comment on the 
important information that neutrino factories could provide. 



1. Introduction 

Partem distribution functions (PDF) are needed to compute hard processes in hadronic 
and nuclear collisions. The method to obtain the PDF from experimental data is well 
established in the case of the free proton: the initial distributions at Qq are evolved 
by the DGLAP equations |Q to larger Q 2 and fitted to available data. The data 
from deep inelastic lepton-proton scattering (DIS) are of main importance in these 
analyses. The nuclear structure functions measured in DIS experiments differ 
from those of the free nuclcons. Dcfinining the ratio vs. deuterium, Rp 2 (x,Q 2 ) = 
^F^(x,Q 2 )/(^F^ ) (x,Q 2 )), several nuclear effects can be distinguished: shadowing 
{Rp 2 < 1) at small values of x, antishadowing (Rp > 1) at intermediate x and EMC 
effect (Rp 2 < 1 again) and Fermi motion at large x. The nuclear effects in F^ translate 
in nuclear PDF (nPDF) which are modified from the ones of the free proton. The 
goal then is to obtain a set of nPDF following the well established procedure used for 
the free proton. In practice, a set of ratios of the PDF in bound and free protons, 
Rf(x, Q 2 ) = ff(x, Q 2 )/ff(x, Q 2 ) for i = g, uv, dv, u. . ., are extracted for a known 
set of the free proton PDF ff(x, Q 2 ). 



2. EKRS analysis 

The main problem in the nuclear case is the lack of experimental data. The DGLAP 
analysis of EKRS, which lead to the set EKS98 ||, uses several sets of DIS data on F^ 
(see Q for the refs.) and data on the Drell-Yan (DY) process measured in pA collisions 
Some other sets of data could be very helpful in constraining the nuclear effects 
for different parton flavours, e.g. charm production to constrain gluons. However, so 
far they are not included in Q because of large error bars (open charm in pA) or the 
presence of final state nuclear effects (charmonium in pA) . Further constraints which 
are used are momentum and baryon number sum rules. At an initial scale, chosen as 
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Qq = 2.25 GeV 2 , the ratios for valence quarks R A (x, Qq) (same for u A and d A ), sea 
quarks R A (x, Qq) (same for u A , d and s A ) and gluons R A (x, Qq) are obtained in the 
following way: 

• At large values of x (cc>0.3) valence quarks dominate, and the data on Rp 2 fix 
the ratio R v but do not constrain the ratio R A . There are no constraints for the 
nuclear gluons, either, in this region. For consistency of the DGLAP evolution, 
it is assumed that R A « R A , and that a similar EMC effect also exists in R A 
already at Qq. 

• At intermediate values of x (0.04<x<0.3) both DIS and DY data constrain 
the ratios R A and R A . The use of DY data || is essential in order to fix the 
relative strength of the valence and sea quark modifications, as DIS alone cannot 
distinguish between them. The baryon number sum rule imposes also constraints 
to R A . The gluon ratio Rf is constrained at 0.02<x<0.2 by the NMC data on 
the Q 2 dependence of the ratio F^/F^ Q, and by momentum conservation. A 
20% antishadowing is found for gluons at x ~ 0.1. 

• At small values of x (x<0.04), F2 is dominated by sea quarks, so DIS data 
constrains mainly R A . The ratio R A is fixed by baryon number conservation 
and turns out to be larger (less shadowing) than R A . At x<0.005, where no 
information from data is obtained in the region Q 2 >1 GeV 2 , a saturation of the 
shadowing (Rp 2 — ► const.) is assumed. This phenomenon has been observed but 
only at Q 2 <c 1 GeV 2 . At the initial scale, R A w Rp 2 is assumed at x<0.01, 
which leads to positive logQ 2 slopes for F A (observed at x ~ 0.01 Q). 

For a given initial condition, LO DGLAP evolution is done. Then, comparing with 
the data at different Q 2 , the best initial distributions are obtained through a recursive 
procedure. The resulting initial ratios at Qq can be seen in Figure [l]. 

3. Comparison with other approaches 

For the moment there is only one global DGLAP analysis on the nPDF similar to 
EKRS §, that of HKM §, §. Figure 1 shows a comparison of the EKS98 and HKM 
results for the ratios R A , R A and R A at Q 2 = 2.25 GeV 2 . The difference between 
the results follows from the fact that the data on the DY process || and on the Q 2 
dependence of F A Q are not used as contraints in the HKM analysis. 

The DY data set is important in the EKRS analysis in fixing R A and R A at 
intermediate x: the DIS data forces R A > 1 at x ~ 0.1 and, as the the DY cross 
sections show almost no nuclear effects (in X2) there, the ratio R A is bound to be less 
than one (no antishadowing for sea quarks). In this conference, preliminary results 
from the HKM analysis with the DY data included were presented 0. As a result, a 
better agreement with EKS98 was found. 

The Q 2 -dependence of the structure function F2 is very sensitive to the gluon 
distribution at small values of x. LO DGLAP evolution of the nPDF gives Q 
dR^, 2 (x,Q 2 )/d\ogQ 2 - R A {2x,Q 2 ) - R^ 2 (x,Q 2 ). NMC has measured g positive 
logQ 2 slopes for the ratio F^/F^- This implies that within the DGLAP framework 
gluon shadowing cannot be much stronger than that in F A in the measured region 
x>0.01. Also too weak a gluon shadowing is outruled || by the NMC data. 
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Valence Sea Gluons 




Figure 1. EKS98 (solid lines) and HKM (dashed lines) nuclear modifications for 
valence, sea and gluon distributions in Pb and C at Q 2 = 2.25 GeV 2 . 



4. Improvements from neutrino DIS data 

DIS experiments with neutrino and antineutrino projectiles could measure M 



F% p = 2x(u + d + s + c) Fl v = 2x (u + d + s + c) 

xF^ p = 2x (— u + d + s — c) X ^3 P = 2x(u — d — s + cj 

with similar relations for neutrons. Different flavors could then be disentangled and 
some of the uncertainties discussed above (e.g. valence at small x and sea at large x) 
would become more directly constrained by data. This would allow for a more detailed 
analysis of the nPDF. Moreover, the valence/sea separation at medium x would be 
measured and could be compared with the results from DY data. This would test 
the universality of the nPDF. Measuring sea and valence quark distributions in vA 
experiments would also shed more light on some open questions of QCD in nuclei, 
such as the probability interpretation of the (n)PDF H. 
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